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Abstract 

In the Eliashberg integral equations for d-wave superconductivity, two different functions 
(a 2 F) n (ui, 9) and (ct 2 F) p ^{uj) determine, respectively, the "normal" self-energy and the "pairing" 
self-energy, oj is the frequency of fluctuations scattering the fermions whose momentum is near 
the Fermi-surface and makes an angle 9 to a chosen axis. We present a quantitative analysis of 
the high-resolution laser based ARPES data on a slightly under doped cuprate compound Bi-2212 
and use the Eliashberg equations to deduce the u and 9 dependence of (a 2 F) n (uj,9) for T just 
above T c and below T c . Besides its detailed ui dependence, we find the remarkable result that this 
function is nearly independent of 9 between the (tt, 7r)-direction and 25 degrees from it, except for 
the dependence of the cut-off energy on 9. Assuming that the same fluctuations determine both 
the normal and the pairing self-energy, we ask what theories give the function (a 2 F) p ^{lo) required 
for the d-wave pairing instability at high temperatures as well as the deduced (a 2 F) n (9, oS). We 
show that the deduced (a 2 F) n (9,uj) can only be obtained from Antiferromagnetic (AFM) fluctu- 
ations if their correlation length is smaller than a lattice constant. Using (a 2 F) p>c i(to) consistent 
with such a correlation length and the symmetry of matrix-elements scattering fermions off AFM 
fluctuations, we calculate T c an show that AFM fluctuations are excluded as the pairing mecha- 
nism for d-wave superconductivity in cuprates. We also consider the quantum-critical fluctuations 
derived microscopically as the fluctuations of the observed loop-current order discovered in the 
under-doped cuprates, and which lead to the marginal fermi-liquid properties in the normal state. 
We show that their frequency dependence and the momentum dependence of their matrix-elements 
to scatter fermions are consistent with the 9 and oj dependence of the deduced {a 2 F) n (tjj,9). The 
pairing kernel (a 2 F) p ^{oj) calculated using the experimental values in the Eliashberg equation 
gives d — wave instability at T c comparable to the experiments. 
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I. INTRODUCTION 



The family of Cuprates have the highest superconducting transition temperatures T c 
discovered so far. Superconductive pairing is in the "d-wave" symmetry. Superconductivity 
with d-wave symmetry and at such high temperatures requires a mechanism different from 
effective electron-electron attraction through virtual exchange of phonons. 

We begin by reviewing the mechanisms for d-wave pairing by a general symmetry analysis 
of the momentum dependence of the pairing vertex. We also define two quantities which 
appear in the analysis of experimental data in the second part and which determine the 
"normal" and the "pairing" self-energies. We then use experimental data from high res- 
olution laser ARPES at different angles across the fermi-surface to deduce aspects of the 
momentum and frequency dependence of the "normal" self-energy both in the normal and 
the superconducting states. We summarize the results of the inversion of the Eliashberg 
Equations to deduce the momentum and frequency dependence of the particle-hole fluc- 
tuations which lead to the observed self-energy. There have been innumerable ideas and 
calculations proposed to understand the properties of the Cuprates. We consider two spe- 
cific mechanisms which have been proposed for d-wave pairing due to a purely electronic 
mechanism: Exchange of Antiferromagnetic fluctuations which are prominent in very un- 
derdoped cuprates and of quantum-critical fluctuations from loop current order observed in 
under doped cuprates in the pseudo-gap region of the phase diagram. We ask, given the 
experimentally deduced fluctuations and their coupling to fermions, whether one or both is 
consistent with the measured ARPES spectrum and the measured T c . 

Soon after BCS theory, Eliashberg pQ used the field theory methodology developed for 
superconductivity by Gorkov [2J to formulate the theory of superconductivity to include 
the frequency dependence of the effective interactions through exchange of phonons. The 
unambiguous experimental proof that the superconductivity in metals such as Pb, Sn, etc. 
is induced by electron-phonon interaction is given by the analysis of tunneling spectrum by 
Rowell and McMillan [3l H] in these metals using the Eliashberg theory and the measurement 
of the spectrum of phonons by neutron scattering. The theory also provides experimental 
proofs of its limit of validity. With one important modification, which does not affect the 
linearized Eliashberg theory which is enough to determine T c , the theory can be used for 
pairing in any symmetry of degenerate fermions due to exchange of any kind of fluctuations, 
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provided the conditions for limit of its validity are satisfied. 

In practice, the procedure for extracting and using information for d-wave superconduc- 
tors and metals, such as the cuprates, in which the high frequency cut-off of the energy of 
relevant fluctuations is an order of magnitude larger than that of lattice vibrations is much 
more demanding of both data and analysis. Even for phonons, Rowell and McMillan used 
data which had a relative accuracy in the measured conductance of 0.2% over a few times 
the highest phonon frequency and from above T c to well below T c . This led to completely 
reliable conclusions. For reasons that have been discussed [5], ordinary tunneling or STM 
is not suitable for deducing spectrum of fluctuations promoting anisotropic pairing as in 
the cuprates. For that we must turn to ARPES and measurements and analysis at various 
angles across the Fermi-surface. The data needs to be consistent and reliable from above T c , 
to at least half of T c , at various angles and over an energy range of about 0.5 eV, which we 
will show is the upper cut-off of the fluctuations. We have relied on the best available laser 
-ARPES data, from the group of one of us (Xingjiang Zhou at Beijing). Even this data at 
present is only reliable in the superconducting state to about 1% up to only an energy of 
about 0.2 eV. The data above T c at angles from the diagonal to the BZ to about 25 degrees 
to it is reliable to this accuracy to about 0.5 eV. We expect future data to completely solve 
this problem but on the basis of existing data some fairly reliable conclusions can be drawn. 

II. PAIRING SYMMETRY 

We will show in this section that in the spin-singlet channel, s-wave pairing is favored 
when the scattering of fermions from k to k' with both near the fermi-surface is independent 
of the angle between k and k' and d-wave pairing is favored when the strongest scattering 
is between k and k' oriented tt/2 with respect to each other, ir/2 scattering can occur both 
through the well known case of the exchange of Antiferromagnetic fluctuations, provided 
they are sharply enough peaked near a commensurate wave-vector which spans the fermi- 
surface, or through exchange of current fluctuations which we will specify below. 

These facts about the favored pairing symmetry are implicit in earlier work jU] and can (in 
most cases) be deduced from the momentum and spin-dependence of the effective interaction 
Hamiltonian, written in terms of the irreducible interaction function -Zs(k, k + q, u) with 
momenta k, k + q at the fermi-surface, total spin S and energy transfer u — > 0. 
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The pairing vertex Is{k,k + q,u) scattering fermions at (k, a), (— k, /3) to 
(k + q, 7), (— k + q, S) with an energy transfer u may be written as 

I s (k, k + q,u) = flf a ,/j(k, k + q)flf 7(i (-k, -k - q)J r a , / 3, 7)( 5(q, k, w). (1) 

^,^,7,15 k, w) is the propagator of the fluctuations which are exchanged by the fermions, 
and <7 aj/ g(k, k + q) is the scattering matrix. 

An interaction of the spin-rotational invariant form may be separated into spin- 
independent and spin-dependent parts: 

I a)/ 3 i7i(5 (k, k') = 1/2 f/i(k,k')5 a ^5 7 5 + J 2 (k, k')cr a/3 • ^^c^c^c-k^Ck^. (2) 

kk' 

The interaction in the spin singlet (S = 0) channel must be even under the interchange 
k — 7- k' and must be odd under this interchange for the triplet (S = 1) channel. Therefore 
both I\ and J 2 can contribute to the S = channel, but only J 2 can contribute to the 
(S = 1) channel. The matrix element for pairing in the S = and S = 1 channels after the 
appropriate spin-traces are 

J(S = 0,k,k') = ^[/ 1 (k,k / )-3(/ 2 (k,k') + / 2 (-k,-k / ))], (3) 
J(S = l,k,k') = i[/ 2 (k,k')-/ 2 (-k,-k')]. (4) 

Consider short-range interactions in real space so that they can be written as separable 
functions of k and k'. The separable functions decompose into sums over different angular 
momentums I if the Fermi-surface is isotropic or more generally into irreducible represen- 
tations of the point-group of the lattice. For pairing on a single Fermi- surface, the physics 
to discern the symmetry of superconductivity can be learnt from considering an isotropic 
fermi-surface with both k and k' on the fermi-surface. In this case case both l\ and J 2 have 
the form ex ^2 e a£(k F )Pi(k)Pi(k'). 

The projected pairing interaction in the £-th angular momentum channel for S = from 
Eq. (J2j, [3]) is usefully written as an integral over the momentum transfer q = |k — k'| ~ 
2kp sin(#/2), where cosd = kp • k' F . Then with sin 8/2 = x and for u — > 0, 

I(S = 0,£) = 2[ dxxP e (l - 2x 2 )[h(2k F x) - 3J 2 (2A; F x)], t even. (5) 
Jo 

Several simple and important points may also be noted from Eq. ([5]): 

(i) An interaction independent of x, i.e a momentum independent interaction gives for all 
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£, except i = 0. This represents the fact that such an interaction is a 5- function in real 
space and any finite I pair wave-function has zero amplitude at the origin. 

(ii) An attractive interaction Ii < is required for the t = case but 1\ < provides a 
repulsive interaction for I ^ 0. 

(iii) An antiferromagnetic I2 interaction is repulsive for i = 0. It is attractive for i = 2 only 
if it peaks sufficiently near the zone boundary. 

(iv) For the contribution both from I\ and from I 2 , I = 2 is favored when the strongest 
interaction is at x 2 = ±1/2, i.e. that the initial and final states are at tt/2 with respect to 
each other. 

It is straightforward to extend this analysis [6] to the more realistic case for the momentum 
dependence of the interaction taking the crystal symmetry into account. We refer to the 
original paper for several different crystal symmetries. For a square lattice, of relevance to 
the cuprates, the peaking of AFM fluctuations near the (n, 7r) point and doping not too far 
from half-filling, strongest scattering of fermions near the Fermi-surface occurs through an 
angle ir/2 leading, as above, to d-wave pairing. At the same time, the "normal" self-energy 
of the fermions depends on their direction of the antiferromagnetic vector Q with respect 
to the crystalline axes. The strength of both effects depends on how steep is the increase 
in I(q) near the zone-boundary and on the details of the electronic dispersion e(k) near the 
chemical potential. 

While the arguments above provide the basic physical principles for discussing pairing 
symmetry, it is necessary also to consider the frequency dependence of the pairing vertex for 
a theory of T c . The considerations determining T c have been discussed elsewhere [7j; here 
we will decipher the spectrum directly from the experiments. 

In order to obtain quantitative conclusions from analysis of ARPES data, it is useful to 
define two functions. From Eq. ([!]) and the ensuing discussion, we define a function which 
enters as the Kernel for Eliashberg Equation for the "pairing" or "anomalous" self-energy 
in the d-wave channel: 



(a 2 F) p4 (u) = [ f ^|^ l( k,k')^(-k,-kO/m^(k,k',c,) (6) 

J \v{k F )\ J \v(k' F )\ v 

- 3^(k,k')^(-k,-k')/mJ- 2 (k,k',u;) N )P 2 (k)P 2 (k')/ f 



Here gi, T\ and #2, J~2 are the spin-independent and spin-dependent vertices and Fluctuation 
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propagators, in complete analogy with Eq. (|2j). P2(k) is the appropriate projection operator 
for the given crystal symmetry, ((l/v2) cos(2#) for isotropic case), and the integrations are 
over the Fermi- surface. 

Equally necessary to define is the Kernel for the equation for the "normal" self-energy in 
which the vertices and the spin-sums enter differently and which explicitly depends on the 
direction on the Fermi-surface: 

J \vQs! F )\ v / 

Note the absolute magnitude signs in ([7]) compared to (pi), because the "normal" self- 
energy must have a pair of conjugate vertices, and the +3, rather than —3 due to spin-trace 
in the particle-hole channel. For s-wave superconductivity, the same function is the kernel 
for both the normal and the pairing self-energies. 

In approximate solutions of the Eliashberg Equations for d-wave pairing, one can define 

K . j^L (8) 

- /^MH (9) 

Here 

(a 2 F(co,9)) n = N(9)(a 2 F) n (k F ,co), (10) 

where N(8) = (\v F (6)\)~ 1 d]t t / 'd6 and k t is the tangent to the fermi- surface. In terms of 
A n , Xd, an approximate solution of the Eliashberg equation yields 

/ l + An \ 

T c ^w c exp l x p^\ (11) 



Better calculations can only be done by numerical solutions but Eq.(|8j) sufficiently explain 
one of the principle points of this paper. 



III. DEDUCED RESULTS FROM ARPES EXPERIMENTS 

In three recent papers [5HTU]. our collaborators and we have presented a fit to the Laser 
ARPES data on a slightly underdoped sample of BISCCO 2212 with T c = 89K and a 
pseudogap temperature of about 125A^. We refer to these papers for the detailed procedure 
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FIG. 1: The Fermi surface of Bi2212 in the first Brillouin zone. The blue solid curve is a calculated 
Fermi-surface and the solid dots are experimentally determined FS at 9 = 0, 5, 10, 15, 20, and 25 
degrees. k± is the distance from the (jr, tt) point. The thick bars along each cut indicate the ranges 
of experimentally measured ARPES MDC data. 

for inverting ARPES ( also given in the nodal direction in Ref. (jTT) in Cuprates. Here we 
only provide a summary of the results. 

Data was available along the momentum cuts shown in Fig. ([I]). We were able to deduce 
the momentum and frequency dependence of the normal state self-energy at T = 107K with 
high accuracy and consistency. The results are shown in Fig. ([2]). The self-energy is labelled 
by the angle 9 with respect to the diagonal of the BZ and by the energy u. No perceptible 
dependence on |k — kp| could be found, except through the upper cut-off in frequency of 
the spectrum, as described below. 

A. Normal State 

Using the deduced self-energy, we have inverted the Eliashberg equation for the normal 
state self-energy to deduce the generalized Eliashberg function (a 2 F(6,e)) n defined in Eq. 
0. The results are shown in Fig. ([3]). 

The deduced (a 2 F(9,e)) n is quite striking in two respects: 

(1) Even though the self-energy, Fig. ^ shows significant angle dependence, the deduced 
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FIG. 2: The real ReT, and imaginary ImTi part of self-energy at T = 107 K for the tilt angles 9 = 
0, 5, 10, 15, 20, and 25 degrees with respect to the diagonal to the BZ and as a function of positive 
energy u. Plot (a) is the result from a detailed fit to band-structure over the whole range of energy 
while (b) is given to compare the results if a linear extrapolation of the band-structure from that 
near the fermi-energy is adopted. 

(a 2 F(9, e)) n is independent of angle to an accuracy of about 10% below an energy of about 
0.2 eV. Above this energy there is an angle dependent cut-off u c (9) for 9 = 20 and 30 degrees 
but the cut-off is independent of angle once it reaches its maximum value of about 0.4 eV, 
as for 9 = 15, 10 and degrees. As has been shown ([§]), the cut-off is at the bottom of the 
band with respect to the fermi-energy at any angle except when this value increases beyond 
about 0.4 eV, where it stays at 0.4 eV. It follows that the spectrum of the fluctuations has 
an intrinsic cut-off which is about 0.4 eV. 

It also follows from the general relation between the self-energy and this fluctuations 
spectra and the bare single-particle Green's function that the angle dependence of the self- 
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FIG. 3: The Eliashberg function (a 2 F) n (8,Lo) discussed in the text extracted from the self-energy 
shown in Fig.Q. Notice that despite the considerable dependence of the self-energy on angle 8, 
the deduced Eliashberg function collapses within the accuracy of its deduction to as single curve 
below about 0.2eV. Above this energy it cuts off, with the cut-off depending on the angle. This 
cut-off is the same, about 0.4 eV for 8 = 0,5, 10 degrees where the bottom of the band is below 
the fermi-energy by larger than 0.4 eV. For 8 = 15, 20, 25 degrees, the cut-off is simply given by 
the measured bottom of the band with respect to the fermi-energy at such angles. 

energy is almost completely due to the dispersion of the bare band. 

(ii) As a function of e, (a 2 F(9, e)) n may be considered to be a sum of two features, a 
nearly ^-independent bump centered at about 0.05 eV and the nearly constant part whose 
intrinsic cut-off, as discussed, is about 0.4 eV. With that value of the cut-off, the 0.05 feature 
has about 10% of the total spectral weight, increasing to about 20% at a cut-off of 0.2 eV. 

The deductions discussed above are consistent with the normal state self-energy deduced 
at 9 = degrees for all the Cuprates by various different groups for which ARPES results 
have been obtained. See fig. Q. The imaginary part of the self-energy for all of them is 
linear in energy (above an energy which depends on temperature or pseudo gap energy) and 
has a distinct energy above which it is constant. This is consistent with scattering from a 
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nearly constant in energy spectrum with a well-defined cut-off. The deduced (a 2 F(9, e)) n at 
9 = is quantitatively consistent with an earlier deduction [11 J from ARPES spectrum in 
the same compound, and also qualitatively consistent with the deduction from optical con- 
ductivity spectrum p2], p3] which preferentially weights the nodal quasi-particles because 
of their larger fermi-velocity. 




Energy (eV) 

FIG. 4: The linewidths of the Momentum Distribution curves for different cuprates as a function 
of the energy near optimal doping and along the diagonal to the BZ. The imaginary part of the 
self-energy is obtained by multiplying this linewidth with the bare fermi-velocity. The detailed 
references for each cuprate are given in Ref. (ll4p 

B. Superconducting State 

As explained above and elsewhere [SHU [9], unlike the Eliashberg Equations for "s-wave" 
pairing, it is necessary in the superconducting state for "d-wave" pairing to extract separately 
the Eliashberg function (a 2 F) n , which occurs as the kernel in the integral equation for the 
"normal" self-energy, and the Eliashberg function (a 2 F) a which occurs as the kernel in the 
integral equation for the " anomalous" or pairing self-energy. The " normal" self-energy has to 
maintain the same sign around the Fermi-surface and therefore is proportional in an isotropic 
system to the form ao + a^ cos 2 (28) + .., or equivalent harmonics in a lattice symmetry, while 
the leading terms for the anomalous self-energy is proportional to b cos(26 l ) +b 2 cos(66 l ) + .... 
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To our knowledge, this separation can only be accomplished only from analysis of Angle- 
Resolved Photoemission data. 

The information on the pairing self-energy is contained only in the difference in the 
ARPES spectra in the superconducting state and the normal state. This difference is ex- 
pected [9] to be less than 1% above an energy of a few times the superconducting gap. Above 
such energy, the noise in the data is at present significantly larger than 1%. Therefore, we 
have not been able to extract the pairing self-energy and to directly deduce (a 2 F) p ^{uj) from 
the data over the fully energy range of 0.4 eV. Below about 0.1 eV, where we can deduce it 
quite well, it essentially follows the angle-independent deduced (a 2 F) n (6, lo). Indeed, we do 
not know of any proposed or reasonable model in which the two functions have quite differ- 
ent frequency dependence in the energy range well above the superconducting gap. But it 
would be ideal to have data which is about 1/2 an order of magnitude better to completely 
settle the shape of (a 2 F(9, e)) a . 

The deduced (a 2 F) n (9, lo) in the superconducting state is noisier than its deduction in 
the normal state, especially for larger 9 and higher energy. Please see ReffJHJ) for the results 
at various temperatures and angles. 

The principal conclusion are that at 9 = degrees, the fluctuations below T c , within the 
uncertainty of determination, are identical to the fluctuations above T c . For larger angles 
and for energies above about 0.1 eV, they are similarly identical to the fluctuations above 
T c . There is a growth of the intensity in the hump around 50 meV for lower temperatures 
and larger angles, and a possible growth of a new peak at about 15 meV. These features 
are probably related to the loss of dissipation due to the opening of the superconducting 
gap and would be interesting to study theoretically in greater detail. However, since they 
cannot affect T c , we will not dwell on them further here. So, in effect we are led only to 
discuss the bulk of the spectra which does not change from below T c to above T c . 

The limit of validity of the Eliashberg theory for the range of parameters we have uncov- 
ered in the Cuprates has been analytically discussed in Ref.([H]). But the conclusive answer 
can only come from analysis of further experiments as we will discuss near the end of the 
next section. 
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IV. DISCUSSION OF RESULTS AND CONCLUSIONS 



In this section, we discuss the implications of the finding that the Eliashberg function 
(a 2 F(6, u)) n defined in Eq. ^ is nearly independent of 6 in the normal state and just below 
T 

- 1 - ti- 
lt should be recognized a priori that the superconducting T c is the property of the normal 
states. Indeed the normal state self-energy just above T c includes all forms of scattering, 
both spin-dependent as well as spin-independent and scattering from all initial momenta k to 
all final momenta k'. It is unlikely that fluctuations invisible in the normal state determine 
T c ; at least no such idea has been expounded. Given this the crucial issue to address is 
how fluctuations which lead to a nearly ^-independent (a 2 F(8,uj)) n are reconcilable with 
the same fluctuations promoting <i-wave pairing, i.e. also giving large values of (a 2 F(u)) aj d, 
defined by Eq. 

We discuss two specific ideas for the source of d-wave superconductivity, (i) Coupling of 
AFM fluctuations to fermions [T3], [TU], [H] , [IE] and (ii) coupling of quantum-critical loop 
current fluctuations to fermions. 

A. AFM fluctuations 

In this section, we show that given the deduced (a 2 F(9,u)) n , AFM fluctuations can be 
excluded as a source of pairing interactions. An important point to note for this and the 
ensuing discussion is that it follows from Eqs. (|6j [7|) that the pairing vertex I depends both 
on the matrix elements g a)i g(k, k + q)<? 7l 5(— k, — k — q) and the spectra of the fluctuations 
Fa,p,'f,6;i'(q.i'k;U n ). The idea of £ = 2 pairing through spin-interactions I 2 [IS], [EZ] 
is quite well developed. For the Hubbard model, the matrix element for scattering the 
fluctuations, p aiJ g(k, k + q}g lt s(— k, — k — q) is momentum independent, being simply the 
square of the local interaction energy U 2 [El [T7j in the weak-coupling limit. At strong 
coupling the matrix element can acquire in addition the symmetry of the exchange energy. 
In neither case can it provide the dominant tt /2 scattering discussed above required for 
d-wave pairing. On the other hand, for a square lattice with AFM fluctuations strongly 
enough peaked at nearly commensurate q, such a scattering is provided by T as discussed 
in Sec. II and elsewhere p2] . 
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Therefore if the fluctuations scattering the fermions are dominantly AFM fluctuation, 
we need to deduce J~ a ,p,j,S;v{^ k, u n ) of an AFM form consistent with the experimentally 
deduced (a 2 F(8, u)) n . From this procedure we can deduce the parameters in the commonly 
assumed form of antiferromagnetic (AF) fluctuations. A phenomenological form for them 
may be written as [HI [19] 

7m^ P (k,k ,«,) = ((k _ k ,_ Q) , e + 1)2 + ( ^ F)2 - (12) 
tuq fixes the integrated spectral weight of the fluctuations, £ is the correlation length and 



uaf, the damping rate of the fluctuations. (a 2 F(6, u)) n can be obtained from Eq.(12) after 
integrating over k. We have followed this procedure to to fit the experimentally deduced 
curve as well as possible, fixing the overall magnitude with m 2 and adjusting the shape with 
the other two parameters. Results of £/a = 1 and 1/ti are shown in Fig. ([5]). While the 
latter gives nearly momentum independent results, the former does so at higher energies but 
not at the lower energies. The discrepancy worsens on increasing £/a or increasing uaf, as 



shown in Fig. (|6j). We then se Eq.(12) with these parameters in the Eliashberg Equations 
(not their approximate solutions in terms of Eq. (|8|) to determine T c . For £/a = 
we find it to be below O.IK for all choices of Uafmi f° r £/ a = 1> T c is calculated as 12 
K for uaf = 0.03eV and 27 K for uaf = O.leV. Only when £/a is larger than 2, do we 
calculate values of T c of order 10 2 K. The experimentally deduced (a 2 F(9, u)) n does not 
allow £/a more than about 1/2. Although no neutron scattering results in the normal state 
are available for Bi-2212 in the normal state, a correlation length £/a < 1 of less than a 
lattice constant is consistent with neutron scattering results in YBa2Cu^0^j r s with 5 near 
the highest T c [20], [2]. 

The conclusion from the above analysis is that if the fluctuations coupling to fermions 
revealed in ARPES experiments were AFM fluctuations, T c would have been less than 10 K. 
The reason for this is clear from looking at Fig (pD. If the correlation length of fluctuations 
is not much larger than the lattice constant, the projection to d-wave scattering, i.e. the 
calculation of the relative dominance of scattering of fermions through 7r/2 is small compared 
to the average or s-wave scattering. Hence of Eq. ([§]) calculated as a (weighted) integral 
over the red curve in Fig. ^ is much smaller than A n which is calculated from the black 
curve. The same conclusions apply to fluctuation spectra of any other operator which is 



of the form of (12) because the data forces it to have a short correlation length. More 
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FIG. 5: The functions discussed in the text for ujaf = 0.03eF. The left panel is for correlation 
length £/a = 1 and the right for £/a = 1/ir. The dashed curves are for from top to bottom from 
the antipodal direction to the nodal directions by increase of 5 degrees. For £/a = l/V, the curves 
are nearly independent of 6. The red curve which determines is a 2 F Pi( [(uj) and the blue curve 
is the s-wave average that determines A n . The T c calculated for £/a = 1/tt and £/a = 1 from the 
solution of the Eliashberg equations are less than 0.1 K and 12 K respectively. 




FIG. 6: The left panel in green shows the result for 0's as in Fig. ^ for £/a = 1 and ojaf = O.leV. 
The red and the blue curves are again as defined in the caption for Fig. §5§. The T c calculated 
with these parameters for AFM fluctuations is 23 K. The right panel shows the variations for the 
same directions for £/a = 2 and ojaf =0.03 eV. The T c calculated for these parameters is 64 K. 



sophisticated correlation functions than (12) cannot change this conclusion significantly. 



B. Quantum-Critical spectra of Loop-Current Order 

We show in this section that a nearly 6 independent (a 2 F(9, u)) n as well as the (a 2 F(uj)) Pt d 
to obtain high T c is consistent with the spectral function of the local quantum-critical fluc- 
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tuations derived recently [22Tf24"] due to the quantum melting of the loop-current order [25J, 
observed universally [2"TH3"U] in the underdoped region of the cuprates, and the coupling 
function of these fluctuations to fermions (31). We briefly summarize the results. 

There are two aspects to the solution: the momentum independence of the derived fluc- 
tuation spectra J-"(q, u) and the momentum dependence of the coupling function of such 
fluctuations to the fermions, g(k, k + q). 

The loop order is specified by a polar time-reversal odd vector L which specifies a pair 
of magnetic fluxes in each unit-cell generated through many body interactions. There are 
four possible orientations of L. In the quantum-fluctuation regime, the orientations of 
L(i,£) in the cells i may be taken to be a continuum. The fluctuations are specified ([22J, 
[23] . [24J by the angular momentum operator U(i, t) which quantum-mechanically flips local 
configurations of the loop order. The Fourier transform of < U + (i, t)U(j, t) > in this theory 
are the propagator of the relevant fluctuations J-(q,u) of Eq.Q. The spectral weight is 
locally critical and of the form, 

-Xotanh(o;/2T), \u\ < uj c ; 
ImF(q,u) = < (13) 

[o, \u)\ >u c . 

This is precisely of the form which leads to the marginal fermi-liquid properties [32] in a 
part of the phase diagram. In the pseudogap region where loop order is observed through 
condensation of the low energy part of these fluctuations, this spectrum is depleted at low 
energies with some of the spectral weight going to weakly momentum dependent collective 
modes ([? ]). Such modes have recently been observed [33] . 

In the continuum limit, U(r) is the angular momentum operator generating rotations 
among the configurations of loop order. Therefore it can only couple to the local angular 
momentum of fermions. So the coupling is of the form 

H mt oc I drj^go ^ + (r,a)(f x p)^(r, <r)U(r) + H.C. (14) 

go is the coupling energy which can be estimated from experiments on single particle scatter- 
ing rate. Hmt has also been explicitly derived for the fermions in a two-dimensional lattice 



model of the Cuprate lattice [51] . It is instructive to note that Eq. (14) is the orbital angular 
momentum analog of the familiar collective spin-fluctuation coupling to spin-flip excitations 
of fermions. 
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We may write Eq. ( 14 ) in momentum space 



Hint =J^9o i(k x k')^ + (k, a)^ + (k', a) ■ U(k - k') + B.C. (15) 

k,k',cr 

This identifies that in the continuum limit (similar conclusions are obtained from the 



lattice version of the calculations), the coupling function g(k, k') of Eq. (14) is given by 



5 (k,k')= 5ol (kxk'), (16) 

giving immediately that the strongest scattering by the fluctuations rotates the momentum 
of the fermions by tt/2. As discussed in Sec. II this promotes "d-wave" pairing if the effective 
interaction is attractive even though the the spectrum of fluctuations itself is momentum 
independent. That it is attractive may be seen as follows: Integrating over the fluctuations 



in Eq.( 15 (gives an effective vertex for scattering of Cooper-pairs: 



A (k, k, ) oc -(k x k') 2 R&F(k - k', u). (17) 

Since ReJ-"(k — k'),uj) < for — uj c < u < u c , independent of momentum, the pairing 
symmetry is given simply by expressing (k x k') 2 in separable form : 

(k x k') 2 = 1/2 [(k 2 x + k 2 y )(k' 2 + k' y 2 ) - (k 2 x - k 2 y )(k' 2 - k'y) 

- (2k x k y )(2k' x k' y )} . (18) 



Pairing interaction in the s-wave channel is repulsive, that in the two d-wave channels is 
equally attractive, and in the odd-parity channels is zero. The factor i in g(k, k'), present 
because the coupling is to fluctuations of time-reversal odd operators, is crucial in determin- 
ing the sign of the interactions of the pairing vertex. In the lattice version of the theory |31j, 
the weighting of the angular decomposition by the density of states in different irreducible 
representations favors the (x 2 — y 2 ) pairing over the 2xy pairing. 

The resolution of the puzzle of a nearly 6 independent a 2 F n (6, u) in Fig. ^ with a strong 
pairing in the d-wave channel comes about in the following way. Consider first the part of 
Fig. ([3]) without the bump. This frequency dependence (including the angle dependent cut- 



off) automatically comes from the derived spectra of Eq. ( 13 ) and this gives no 6 dependence 
in a 2 F n (6,uj), which comes only from integrating |g(k, k')| 2 over k' as in Eq. Q and then 
converting the resulting function of k to that of 6 by multiplying with the local density of 
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FIG. 7: The angle dependence for any frequency expected for a 2 F n (9,co) in the coupling of loop- 
order fluctuations to fermions in a simple model for the matrix elements. The red-curve is for the 
ratio of velocity in the (tt, tt) direction to that in the (tt, 0) direction of 2 and the black curve for 
4. The red-curve is slightly displaced with respect to the black curve for clarity. 



states at the chemical potential N(6). In the continuum limit Eq.(16), this automatically 
yields a ^-independent result. We have also calculated it from the leading lattice-dependent 
g(k, k'), given in Ref.j 



g(k, k') oc igo ( sin(k x a) sm(k' y a) — sm(k y a) (k' x a)) (19) 

The results are shown in Fig. ^ for two choices of go, which determines the ratio of the 
fermi-velocity in the (tt, it) directions to that in the (tc, 0) directions to be 2 and 4. The 
experimental result for this ratio is about 2. The variation in the 9 dependence in the 
measured range, tt, tt direction to about tt/8 from it is about 10%, which on comparison 
with Fig. ([3]) is within the experimental uncertainty. With higher resolution data, the 
prediction is that there should be a frequency-independent ^-dependence of the order of 
magnitude given in Fig. . 

There is one aspect of the deduced a 2 F(9, u) which is not given by the critical theory. 
This is the low energy bump as a function of u. If the bump occurs only in samples studied in 
the pseudogap region, it is probably due to the collective modes which have been discovered 
to be special to this region of the phase diagram [531 El] • This can be checked by equally 
high resolution data in samples in the quantum-critical region or the overdoped region of 
the phase diagram. Alternatively, it may be due to coupling to oxygen vibrations modes, as 
has been suggested [35] . 

In the loop-fluctuation theory, the equivalent of the red and the black curves of Fig. (§ 
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used in the Eliashberg Equations to calculate T c have identical frequency dependence and 
nearly the same magnitude. The calculated T c using the experimentally deduced a 2 F(9, u) 
is about 100 K. There remain valid questions about the applicability of the Eliashberg 
equation to spectra whose cut-off is comparable to the over-all electronic bandwidth and 
with coupling constants of about 1/2. We believe that the best way to investigate this is to 
deduce the variations in a 2 F(8, u) for a variety of dopings and T c in a given compound. We 
doubt that Eliashberg theory works quantitatively to the accuracy that it does for electron- 
phonon interaction but we expect to get the ratios of calculated T c 's for different dopings to 
compare well with the ratio of observed T c 's. Such work is in progress. 



V. CONCLUDING REMARKS 




FIG. 8: The figure on the left shows the proposed [36] Universal Phase diagram for hole-doped 
cuprates. The right part of the figure (taken from Ref.([37J) shows the lines of phase transitions 
from experiments in Hgl201(red) and YBCO(blue) at low and intermediate dopings. The main 
panel and the inset show data plotted with two different way of determining p. Stars denote T c , 
and circles denote T* , the temperature below which the predicted magnetic order is seen. 

The phase diagram in Fig. ^ may be taken to provide a motivating principle in the 
search for a theory of the Curpates - that the universal properties of all the regions in it be 
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governed by a single idea which can be calculated in a controlled fashion for properties in 
all the regions and with predictions which are testable in experiments. 

The left part of Fig. ^ divides the T — doping plane of hole-doped cuprates into re- 
gions with different characteristic thermodynamic and transport properties. Such a diagram 
was first proposed in 1997 [36, 38J, but only in the past five years have some of its princi- 
ple features been verified in experiments. The essential feature of the phase diagram is a 
quantum-critical point from which emanate three regions - a region II of broken symmetry, a 
region I whose properties are determined by the quantum-critical fluctuations of a proposed 
order in region II, with a crossover to a fermi-liquid region III. D-wave superconductivity 
straddles the three regions at low temperatures. A time-reversal breaking order, preserving 
translational symmetry was predicted for the region II. Such an order has been observed by 
polarized neutron scattering in YBa 2 Cu 3 6+x [39J , HgBaCuO gO], LaSrCuO [H], BISCCO 
[12], the last consistent with earlier dichroic ARPES experiments [13]. The second part of 
Fig. [8] taken from Ref.([37]) shows the results for T c and the temperature of the predicted 
order, which coincides within uncertainty with T* determined by other methods for the first 
two compounds. In the underdoped region, other orders have been discovered; they are dif- 
ferent in different cuprates and do not demarcate the universal change in properties as the 
lines drawn in Fig. ([8]) do. For example the sharp decrease in the order in Y Ba^Cu-^O^x 
at small x coincides with the onset [H] possibly of smectic or stripe order. 

The motivation for proposing a symmetry breaking in the pseudogap regime was 
the "strange-metal" Region I, whose properties could be understood by hypothesizing a 
quantum-critical spectrum [32J. Two major theoretical developments since then have been 
the microscopic derivation of this quantum critical spectrum as due to the fluctuations of 
the order parameter in Region II, and the derivation of the matrix element for coupling the 
fluctuations to fermions and the proof that it promotes d-wave pairing as discussed above. 

The different aspects of the theory: the loop-current order, its quantum-critical fluctu- 
ations and their coupling to fermions are all essential and used in the theory of d-wave 
superconductivity summarized above. We have shown that the detailed analysis of high 
resolution ARPES yields results which are very specific in their angle and frequency depen- 
dence and that they are consistent for T c in d-wave symmetry from the same theory with 
which aspects of the normal state in Region I are understood and with which the order ob- 
served in Region II was predicted. There are aspects of experiments in Region II, such as the 
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single-particle spectra and the magneto-oscillations which however remain to be understood 
on the basis of this (or any other) ideas. 

We have also shown that any existing theory based on momentum-dependent fluctuation 
spectra such as AFM fluctuations would need to have a correlation length smaller than 
the lattice constant to be consistent with the deductions from ARPES and that such small 
correlation lengths would yield T c in the d-wave channel with too low a value. This is hardly 
surprising since theories based on such ideas have not succeeded in understanding the normal 
state properties in the Cuprates. 

Acknowledgements: Thanks are due to Elihu Abrahams for a critical reading of the 
manuscript. CMV's research is supported by a grant from the Division of Materials Research 
of the National Science Foundation. 



[1] G. M. Eliashberg. Zh. Eksp. Teor. Fiz., 38:966, 1960. 

[2] L. P. Gorkov. JETP, 7:505, 1960. 

[3] W. McMillan and J.R. Rowell. Phys. Rev. Lett, 14:108, 1965. 

[4] J. R. Schrieffer, D. J. Scalapino, and J. W. Wilkins. Phys. Rev. Lett., 10:336, 1963. 

[5] I. Vekhter and C. M. Varma. Phys. Rev. Lett., 90:237003, 2003. 

[6] S. Schmitt-Rink K. Miyake and C. M. Varma. Phys. Rev., B 34:6554, 1986. 

[7] CM. Varma. arXiv.1001.3618. 

[8] J.M. Bok et al. Phys. Rev. B, 81:174516, 2010. 

[9] J.H. Yun et al. Phys. Rev. B, 814:104521, 2011. 

[10] Wenten Zhao et al. Phys. Rev. B, to be published. 

[11] E. Schachinger and J.P. Carbotte. Phys. Rev. B, 77:094525, 2008. 

[12] J. Hwang et al. Phys. Rev. B, 75:144508, 2007. 

[13] E. van Hemmen et al. Phys. Rev. B, 79:184512, 2009. 

[14] Lijun Zhu et al. Phys. Rev. Lett, 100:057001, 2008. 

[15] K. Miyake, S. Schmitt-Rink, and CM. Varma. Phys. Rev. B, 34:6554, 1986. 

[16] D. J. Scalapino, E. Loh, and J. E. Hirsch. Phys. Rev. B, 34:8190, 1986. 

[17] D.J. Scalapino. Physics Reports, 250:329, 1995. 

[18] H. Monien A.J. Millis and D. Pines. Phys. Rev. B, 42:167, 1990. 



21 



[19] D. Pines A.V. Chubukov and J. Schmalian. The Physics of Superconductors- vol. 2. Springer, 
Berlin, 2003. 



[20 
[21 
[22 
[23 
[24 
[25 
[26 

[27; 

[28 
[29 
[30 
[31 
[32 
[33 
[34 
[35 
[36 
[37; 
[38 
[39 
[40 
[41 
[42 
[43 
[44; 



H.A. Mook et al. Phys. Rev. Lett, 70:3490, 1993. 

A. V. Balatasky and P. Bourges. Phys. Rev. Lett, 82:5337, 1990. 
V. Aji and C. M. Varma. Phys. Rev. Lett, 99:067003, 2007. 

V. Aji and C. M. Varma. Phys. Rev. B, 79:184501, 2009. 

V. Aji and C. M. Varma. Phys. Rev. B, 82:174501, 2010. 

M. E. Simon and C. M. Varma. Phys. Rev. Lett, 86:1781, 2001. 

C M. Varma. Phys. Rev. B, 73:155113, 2006. 

B. Fauque et al. Phys. Rev. Lett, 96:197001, 2006. 
H. Mook et al. Phys. Rev. B, 78:020506, 2008. 

Y. Li et al. PNature, 455:372, 2008. 

A. Kaminski et al. Nature, 416:610, 2002. 

A. Shekhter V. Aji and CM. Varma. Phys. Rev. B, 81:064515, 2010. 
CM. Varma et al. Phys. Rev. Lett, 63:1996, 1989. 

Y. Li et al. Nature, 468:283, 2010. 

Yan He and CM. Varma. Phys. Rev. Lett, 106:147001, 2011. 
S. Johnston et al. Phys. Rev. B, 82:064513, 2010. 
CM. Varma. Phys. Rev. B, 55:14554, 1997. 
Y. Li et al. Phys. Rev. B, 2011. 

J.L. Tallon and J.W. Loram. Physica C, 349:53, 2000. 

B. Fauque et al. Phys. Rev. Lett, 96:197001, 2006. 
Y. Li et al. Nature, 455:372, 2008. 

V. Baledent et al. Phys. Rev. Lett, 105:0207004, 2010. 
P. Bourges (private comm.). to be published, 2010. 
A. Kaminski et al. Nature, 416:610, 2002. 

P. Bourges and Y. Sidis. Comptes Rendus Physique, 12:461, 2011. 



22 



